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Abstract 

Mitochondrial DNA (mtDNA) mutations can result in mitochondrial dysfunction, but emerging experimental data question 
the fundamental role of mtDNA mutagenesis in age-associated mitochondrial impairment. The multicopy nature of mtDNA 
renders the impact of a given mtDNA mutation unpredictable. In this study, we compared mtDNA stability and mtRNA 
integrity during normal aging. Seven distinct sites in mouse brain mtDNA and corresponding mtRNA were analyzed. 
Accumulation of mtDNA mutations during aging was highly site-specific. The variation in mutation frequencies overrode 
the age-mediated increase by more than 100-fold and aging generally did not influence mtDNA mutagenesis. Errors 
introduced by mtRNA polymerase were also site-dependent and up to two hundred-fold more frequent than mtDNA 
mutations, and independent of mtDNA mutation frequency. We therefore conclude that mitochondrial transcription fidelity 
limits the impact of mtDNA mutations. 
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Introduction 

The mitochondrial theory of aging postulates that mutations in 
mitochondrial DNA (mtDNA) accumulate with age and result in 
impaired quality and activity of the mtDNA-encoded proteins. 
The theory is supported by the fact that mitochondrial function 
decreases with age, presumably due to accumulation of somatic 
mtDNA mutations [1]. Mitochondrial dysfunction associated with 
accumulation of clonal expansions of deletion mutations are 
reported in nucleoside reverse transcriptase inhibitor (NRTI)- 
treated individuals [2,3] and neurodegenerative disorders includ- 
ing Multiple Sclerosis, Alzheimer's Disease and Parkinson's 
Disease [4]. Although mutated mtDNA are selectively removed 
during foUiculogenesis [5] as well as during maternal transmission 
[6], inheritable heteroplasmy causes inheritable mitochondrial 
disease like MELAS when such processes fail [7]. 

The Polg mutator mouse expresses an error-prone mtDNA 
polymerase f and accumulates excessive mtDNA mutations in an 
age-dependent manner. The Polg mutator evidently demonstrates 
that mtDNA mutations can result in pathology and shortened 
lifespan [8] . Although this model has been used to demonstrate the 
correlation between mtDNA mutations, mitochondrial dysfunc- 
tion and premature aging, it is questionable to what extent this 
genetic mitochondrial mutator model represents the molecular 
mechanisms that underlie mitochondrial dysfunction during 
normal aging. For instance, the level of mtDNA substitution 



mutations in old individuals with a corresponding mitochondrial 
dysfunction does not produce a mitochondrial dysfunction when 
present in young mutator mice [9]. Rather, mtDNA deletions 
were suggested to be responsible for age-mediated dysfunction. 
Accumulation of mtDNA deletions correlated with the phenotype 
and mutation frequency during normal aging [10]. The deleteri- 
ous functional effects of deletion mutations are demonstrated 
clinically in Kearns-Sayre Syndrome [11]. However, genetic 
models have demonstrated that deletion mutations in up to 60% of 
the mtDNA molecules are tolerated without manifesting into 
phenotypic abnormalities [12]. In view of the relatively high 
tolerance for mtDNA deletion mutations, there is an unexplained 
discrepancy between the observed mutation frequency during 
normal aging and the age-associated dysfunction [12-14]. 

The functional impact of mtDNA mutations is hard to predict 
because of the multiplicity of mtDNA molecules in the cell. 
mtDNA copy number is also subjected to variations and the 
heteroplasmic state implies that the likelihood for a mutation to 
manifest into dysfunctional protein depends on the mtDNA copy 
number, given that transcription occurs randomly among mtDNA 
molecules. In view of the redundancy of mtDNA to serve 
templates for downstream mitochondrial protein components, 
we reasoned that the best strategy to evaluate mtDNA mutagenesis 
would be to address the integrity of the mitochondrial RNA 
(mtRNA). In order for mtDNA mutations to result in functional 
impairment, the mutations must significantly modify the popula- 
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tion of mtRNA molecules. To investigate the impact of mtDNA 
mutations with age, we developed an assay to determine mtRNA 
integrity with high resolution and used this technology to compare 
mtDNA mutagenesis and mtRNA error frequency in brains from 
young mice with those from old mice, which were associated with 
impaired mitochondrial function. Our results show that mtRNA 
error frequency can be used to validate mtDNA mutagenesis. 

Materials and Methods 

Ethics statements 

The project was approved by the Norwegian Animal Research 
Authority and in accordance with the laws and regulations 
controlling experimental procedures in Norway and the European 
Union's Directive 86/609/EEC. 

Animal and Tissue Collection 

WUd-type C57BL/6 mice (males /females) were either bred in- 
house (old group, 18 months, n = 8) or purchased from Taconic 
(young group, 1 month, n = 8). Heterozygous (Polg™"'^''^, n = 3) 
and homozygous (Polg™"'^™"', n = 3) mutator mice were about 9 
months old (36-39 weeks). The project was approved by the 
Norwegian Animal Research Authorit}' and in accordant^ with 
the laws and regulations controlling experimental procedures in 
Norway and the European Union's Directive 86/609/EEC. 
Mouse whole brain without cerebellum was collected immediately 
after the animals were sacrificed. One of the hemispheres was 
submerged in RNAlater solution (Sigma) for DNA/RNA extrac- 
tion and the other used for detection of mitochondrial functions. 

Isolation and Quantification of Nucleic Acids 

The RNAlater-stored brain tissues were homogenized with 
FastPrep-24 (MP Biomedical, USA). Total DNA was isolated by 
the DNeasy Blood and tissue kit (Qiagen). Total RNA was isolated 
using the RNeasy Mini Kit (Qiagen) according to manufacturer's 
protocols with modifications. In brief, 20-30 mg of homogenized 
brain tissue was lysed and applied to spin column. After the 
column was washed with AWl buffer, DNase I (Qiagen) was 
added to the column and incubated at room temperature for 
15 min to remove residual DNA. To completely remove possible 
contamination of the residual DNA, eluted RNA was further 
treated with rDNase I at 37 °C for 20 min according to 
manufacturer's protocol (DNase treatment & removal, Ambion). 
Quantification of DNA and RNA was done on Nanodrop 
spectrophotometer (Thermo Scientific). 

Analysis of mtDNA Point Mutation Frequency 

Double stranded mtDNA mutation frequency was estimated as 
described previously with modifications [15]. Briefly, total DNA 
was treated with SI nuclease (Qiagen; 10 U) for 15 min at 37 °C 
and subsequently digested with TaqI restriction enzyme (New 
England Biolabs; 100 U) at 65 °C for 1 h to remove single strands, 
damaged DNA and non-mutated DNA. The remaining mutated 
TaqI restriction sites were quantified in a qPCR reaction using the 
specific primers listed in Table SI in File SI. To ensure complete 
digestion, an additional TaqI treatment (100 U; 65 °C for 1 5 min) 
was performed on the digested DNA samples prior to the 
subsequent qPCR analyses. Complete digestion was verified by 
subsequent TaqI treatment of qPCR products. Mutation frequen- 
cy is found as (2exp(CT^'"'' " CT ^'^) * 4) ' per nt. The presence 
of nuclear mtDNA fragments (numts) [16] could potentially 
interfere with our data as we analyze total DNA and pseudogenes 
with altered nucleotide sequence potentially could generate 
positive signals in our restriction cleavage inhibition-based assay. 



However, the melt-curve analyses indicated that we always 
amplified the same PGR fragment. The CT values for the seven 
mtDNA sites were in the same range and thus supports that do not 
interfere with our results. 

Analysis of mtRNA Errors 

cDNA was prepared from 0.5-1.0 ng total RNA using the High 
Capacity cDNA reverse transcription kit (Applied Biosystems). 
Subsequently, in order to create- double-stranded cDNA for TaqI 
digestion, PCR was performc-d with 100 ng of cDNA with a 
HotStart PCR system (5 Prime) in a total volume of 50 |.ll 
containing 1 x HotStart PCR bufier, 0.2 mM of dNTPs, 0.2 ^M 
of each specific forward and reverse primer and 2.5 U of Taq 
DNA polymerase. The PCR reactions were performed for 10 
cycles. The amplified products were diluted to give appropriate 
initial input for qPCR quantification. Then 100 ng of DNA were 
either left untreated (nt) or treated with TaqI (TaqI; 100 U, 65 °C 
for 1 h). The mtRNA frequency was calculated from the resulting 
CT values using the same formula as for mtDNA mutations. The 
procedure is outhned in Figure SI. 

The mtRNA deletions were evaluated by analyzing possible 
truncated products in a PCR amplification of extended cDNA 
regions. cDNA was first amplified using primers 12S-F6 (Table SI 
in File SI) (full-length product: 685 bp trancripts) of 12S ribosomal 
RNA with a HotStart PCR for 19 cycles in order to obtain 
appropriate amounts of DNA, whereupon an aliquote of 1 00 ng 
was digested with TaqI enzyme (100 U, 65°C for 1 h). 6 ng of 
digested or non-digested samples were further used as templates to 
ampKfy two amplicons of 206 and 442 bp with primers 12S and 
12S-F4 in a 40-cycle PCR reaction, respectively (Table SI in File 
SI). The amplified products were analyzed by Agilent 2100 
Bioanalyzer using Agilent DNA 1000 regents (G2938-80020, 
Agilent Technologies). 

Mitochondrial Gene Expression 

Mitochondrial gene expression analyses were performed by 
qPCR with the primer sets provided in Table SI in File SI. Total 
RNA was isolated and reverse-transcribed into cDNA and 
quantified by qPCR. Expression levels were determined from 
the difference in ACT values. The housekeeping gene Gapdh was 
used as internal control. 

Western Analysis 

Mitochondrial proteins (5-10 |Xg) were separated on a 12% 
NuPAGE Bis-Tris gel (Invitrogen) and blotted onto PVDF 
membranes (Millipore). After blocking in 5% dry milk in TEST 
buffer (50 mM Tris pH 7.4, 150 mM NaCl, 0.1% Tween 20) for 
1 h at room temperature, blots were probed for specific subunits of 
the respiratory complexes overnight at 4°C. Primary antibodies 
used were: anti-total OXPHOS antibody cocktail (1:2000, 
MitoSciences), anti-ND6 (1:400, Biorbyt) and as loading control 
anti-VDAC (1:500, Abeam). As secondary antibodies alkaline 
phosphatase (AP)-conjugated goat anti-rabbit IgG (1:30000) and 
goat anti-mouse IgG (1:5000) from Santa Cruz Biotechnology 
were used. Fluorescent signals from enhanced chemifluorescence 
substrate (GE Hesilthcare) were detected by Bio-Rad ChemiDoc 
XRS+ and analyzed with Image Lab Software (Bio-Rad). 

Statistical Analysis 

Pearson correlation was used to evaluate correlations between 
DNA and RNA and SPSS version 10 was used for the statistical 
analysis. WUcoxon signed-rank test was employed to evaluate 
effect of age on mutagenesis and RNA integrity. Student's t-test 
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Figure 1 . Site- and Age- Dependent Accumulation of Mutations In Brain mtDNA. (A) DNA from young (1 month, n = 8) and old (1 8 months, 
n = 8) mice were analyzed for mutation rate in 7 different sites in the mtDNA. DNA was isolated from whole brain. (B) DNA from heterozygous (mut/+, 
n = 3) and homozygous (mut/mut, n = 3) mutator mice (9 months) were analyzed in the same sites. Figures show mean with SD. **p<0.01 , *p< 0.05. 
doi:1 0.1 371/journal.pone.0096940.g001 



with two-tailed distribution and equal variance was used to 
calculate significance, *p<0.05; **p<0.01. 

Results 

mtDNA Mutation Frequency Varies More Between 
Specific Sites Than With Age 

We used a qPCR-based method for high-resolution quantifica- 
tion of mtDNA mutations in TaqI restriction sites that was 
vahdated and described previously [15]. Total DNA from the 
hemispheres was used for the analyses rather than DNA from an 
isolated subfraction of brain mitochondria. We selected seven 
dilferent sites in the mitochondrial genome for assessment of 
mutations in mtDNA during aging. We chose mice of 1 and 18 
months age to represent young and old animals, respectively. 
Although the age of 1 month does not represent all aspects of a 
developed mouse brain, we selected this age to limit putative 
contribution of age on somatic mutagenesis. 18 months was 
likewise chosen to avoid putative age-associated disease-driven 
mtDNA mutagenesis from influencing the results. We confirmed 
that brain mitochondrial function was impaired in the old cohort 
(about 35% reduction in complex I, II and V; Figure SI). 

The mtDNA mutation frequencies in brains from young 
individuals varied remarkably between the different loci, from 
1.3 X 10"*^ per lit in the 1 2S ribosomal RNA to 188 x 10"'^pernt 
in the Coxl site (Figure lA). Mitochondrial DNA from aged brains 
(18 months) displayed significantly elevated mutagenesis in the 
NdS (+ 40%) and Ndl (+ 90%) sites and borderline significant 
(p = 0.07) increase in the 12S site (+ 50%), whereas the NdS, Nd6, 
Cytb and Coxl sites were unaffected. Mutations did not generally 
increased with age as implied from the WUcoxon signed-rank test. 
In line with the common responsiveness of the Nd 1 and NdS sites 
to age, we observed a significant correlation between these 
(r = 0.58) among the individuals. A significant correlation was 
identified between 12S and NdS (r = 0.61) and 12S and Nd6 
(r = 0.744) as well. These data suggest that increased mtDNA 
mutagenesis during age is region dependent. We included DNA 
from Polg mutator mice that express one or two alleles of the 3 '-5' 
exonuclease-deficieiit DNA polymerase y as controls. The mutator 
mice accumulate several hundredfold more mutations than control 



mice with age [8] . This was confirmed by us as mtDNA from the 
Polg mutator mice (9 months age) contained up to 400-fold higher 
levels of mutations than control mice, dependent on the site and 
allelic dose (Figure IB). There was no apparent correlation 
between the site-specific mutation level and the effect of 5' 
exonuclease-deficiency. However, the 12S, Ndl and NdS sites, 
which were relatively most affected in the Polg mutator mice (40 to 
400-fold versus < 5-fold for the remaining) were those that 
increased significantiy with age. 

Site-Specific Transcriptional Errors are Independent of 
Age and Mutation Frequency 

The large variations in mutation frequency suggested that the 
penetrance might be site-specific. Since there are multiple copies 
of mtDNA, the frequency of mutations must be higher than the 
rate of ribonucleotide misincorporation by mitochondrial RNA 
polymerase at the specific site in order to manifest into 
mitochondrial dysfunction. We developed a high resolution assay 
to analyze mtRNA errors introduced in the same TaqI restriction 
recognition sequences. The procedure is illustrated in Figure 2A. 
Total RNA was isolated from the same samples as used for DNA 
isolations, and used in reverse transcriptase reactions to produce 
cDNA that was subsequently amplified by PCR and finally 
analyzed for Taql-resistant sequences. To control for the artificial 
errors introduced by this procedure, we compared misincorpora- 
tion of two high fidelity reverse transcriptases and estimated 
incorporation errors by the DNA polymerase in the subsequent 
PCR amplification. RNA from the mutator mouse served as an 
ultimate control for RNA errors that manifest from mutations in 
mtDNA. Two commercially available high fidelity reverse 
transcriptases (with a declared fidelity of 2 to 7 x 10 ' per nt 
according to the manufacturer) gave identical error rates in the 
Coxl site for the different genotypes, demonstrating that the 
fidelity of reverse transcriptases is not restricting the conclusions 
(Figure 2B). The 12S and Ndl sites were additionally tested by the 
two reverse transcriptases and gave similar conclusion (data not 
shown). Furthermore, by amplifying DNA and subsequentiy 
quantifying mutated PCR products, we determined that the errors 
introduced by 6 cycles of PCR (from 9 to 15) yielded less than 2 x 



PLOS ONE I www.plosone.org 



3 



May 2014 | Volume 9 | Issue 5 | e96940 



mtRNA Errors Override mtDNA Mutations 



B 



Tissue collection 



i 



Rapid isolation, store ir> RNAlater 



Total RNA isolation 



i 



D3 JO e DN=5; t'eafneit '5qjT5d to 



cONA synthesis random priming 



i 



:DN-. 3viti5S3 by jjc of '5od3n 



Linear PGR amplification 
100 ngcDrjA; 10 cycles 



i 



Stachionrtetrt: amplif cation from t 1 

'n' "'Ol :20 5S of jlS'i friir-^ots 



Quantify PGR products 



i 



qPCR =1= vS3 of sjojsqjiit PCR 
o'od j:t 



Prepare required DNA concentration 
(20 ng/ml;equivalentto0.15nlVI) 



i 



Required to fsciiitae rsprodudbls 
'sst' :t 31 t'ssfneit 



TaqI restriction digestion 
(100ngONA:lhat65=G) 



i 



?5'f3"r! 10='= s \'. ti =id wrtnout 
T=q! 'i;t' :t 4i:v"^i 



qPGR analysis 



i 



Rji t' 0 :=ts2 t3 obtan Sv-'sg 



RNA integrity 



6000 



o 4000 



1- o 
o ^ 



< 

Z 

•i—t 

E 



2000 



c 

o 



c 
o 



I Applied Biosystems 
I Stratagem 




IM 



18M mut/+ mut/mut 



15 -1 



10 



5 - 



4 



RNAintes'-tv = 1 



9 C7C 



12 eye 



15 eye 



Figure 2. Validation of Method. (A) Flow chart illustrating the procedure for detection of mtRNA integrity. (B) Comparison of RNA integrity in Coxl 
site obtained by using two different commercially available reverse transcriptases. (C) Estimation of errors introduced in the PCR amplication. The 12S 
rRNA region was amplified with increasing cycle number and 100 ng PCR product was either digested with TaqI or left untreated and subsequently 
analyzed for mutations. The mutation frequency is plotted as function of PCR cycle number. The figures show mean with SD from three independent 
experiments. 

doi:10.1371/journal.pone.0096940.g002 



10 ^ per nt (Figure 2C). Together, the procedure therefore 
generates artificial errors with an estimated fi-equency of less than 1 
X 10 per nt. When this procedure was applied on samples from 
young and old mice, we found that the frequencies of 
misincorporated ribonucleotides in mtRNA in vivo were at least 
3-fold higher than this and varied in a site-specific manner 
(Figure 3). Overall, there was a 35-fold variation in error 
frequency. Similarly as for mtDNA mutation frequency, we used 
WUcoxon signed-rank test on relative average mtRNA error levels 
and fovmd that the two age groups were not significantly different 
(p>0.05). The mtRNA error frequency spans from 0.3 to 15 x 
10~* per nt, and overrides site-specific mtDNA mutation 
frequency from 3 to 200-fold (Figure 3A). We did not find any 
correlation between mtRNA error frequencies and mtDNA 
mutations, except for an interesting pattern in the Nd5/Nd6 sites 
where the mutation frequency in Nd6 displayed a negative 
correlation with mtRNA integrity in the neighboring Nd5 site 



(r=— 0.59, p = 0.021). The lack of positive correlation between 
mutations and mtRNA errors strongly indicates that the mtDNA 
mutation frequency at a given site does not predict the 
corresponding mtRNA error frequency. Only mtDNA from the 
Polg mutator mice contained mutations at a level that possibly 
could influence the mtRNA integrity (Figure IB). This was verified 
by Pearson correlation analyses, which demonstrated strong 
correlation (r>0.84) between mtDNA mutation and mtRNA 
errors in the 12S, Ndl, Nd3, Cytb and Coxl sites in the mutator 
mice, thereby demonstrating the penetrance of the mutator 
genotype. The homozygous mutator mice had in average 2.6- 
fold more mtRNA errors than wt. In contrast, the mtRNA error 
frequencies in the heterozygous mutator mice were not signifi- 
cantiy different from wt, except for the Ndl site (Figure 3B). The 
relatively low effect of single allele PolG mutation on mtRNA 
errors in heterozygous mutator mice correlates with its innocuous 
phenotype compared to homozygous mutator mouse [8]. 
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Figure 3. RNA Error Frequency is Age-independent in Controls but Elevated in Mutator Mice. (A) mtRNA errors were measured as 
described at the same 7 loci as for mtDNA from brains of young (1 months, n = 8) and old (18 months, n = 8). (B) mtRNA errors in heterozygous (mut/+ 
, n = 3) and homozygous (mut/mut, n = 3) mice. Figures show mean with SD. p**<0.01 vs. 18 months, p*< 0.05 vs. 18 months. 
doi:10.1371/journal.pone.0096940.g003 



Since the relatively low mtRNA error frequency determined for 
the Nd6 site is in the range of the estimated resolution of the 
method, it was a possibility that the Nd6 mtRNA error frequency 
was overestimated. As an internal control of the resolution of the 
assay, we addressed the frequency of errors in nuclear ribosomal 
RNA by the same procedure. However, in contrast to mtRNA, we 
were not able to quantify the exact number of Taql-resistant 
cDNA from the nuclear 28S ribosomal RNA because the amount 
of potential templates was under the detection limit. However, 
based on the amount of non-digested molecules and this limit, we 
estimated an upper range of the errors in the 28S rRNA to be 2 x 
10"'' per nt (Figure S2). The high fidelity of 28S rRNA 
transcription is expected because of the strong proofreading 
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Figure 4. Level of RNA Deletions is Age-Independent. Intentional 
capture of mtRNA deletions by PCR amplification of extended cDNA 
regions. 

doi:1 0.1 371 /journal.pone.0096940.g004 



Figure 5. Stable Expression of Genes and Proteins with Age. (Aj 

Expression of mitochondrial genes was analyzed by RT-qPCR normal- 
ized to Gapdh and presented relative to 1 month old mice. (B) 
Expression levels of proteins were evaluated by western analyses using 
antibodies as specified in Experimental Procedures. The relative 
increase in COXI and ND6 with age is indicated. Figure shows mean 
with SD. **p<0.01, *p< 0.05. 
doi:1 0.1 371/journal.pone.0096940.g005 
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activity of RNA polymerase I [17] Hence, these results show that 
the resolution of the procedure is sufliciently high to detect 
mtRNA errors but not rRNA errors, and demonstrates that the 
lowest error frequency in the Nd6 site is 3.2 x 10~^ per nt. 

No Indications of mtRNA Alterations with Age 

We evaluated the impact of mtDNA deletion mutations by an 
analogue strategy. The contribution of mtDNA deletions to errors 
in mtRNA was determined by intentional capture of mtRNA 
deletions by PGR amplification of extended cDNA regions. We 
did not identify detectable deleted transcripts in the template 
mixture by this approach (Figure 4). After removing intact cDNA 
by TaqI digestion prior to PGR capturing, we observed truncated 
PGR products from both age groups. However, despite the 
preference for truncated products to be amplified in PGR, full- 
length PGR products were still more abundant in both young and 
old brain. We conclude that deletions contribute less than single 
base substitutions to the TaqI resistance (Figure 4). The innocuous 
impact of age on mtRNA integrity was underlined by the normal 
or increased mitochondrial gene expression in the old group 
(Figure 5A), which correlated with elevated expression of 
mitochondrial proteins in the older individuals (Figure 5B). 

Discussion 

Here, we show that subtle increases in mtDNA mutation 
frequency during normal aging are insufficient to influence the 
quality of mtRNA in seven selected sites. The mtRNA integrity is 
determined by the fidelity of the mtRNA polymerase, which masks 
mtDNA mutation frequency by an order that is dependent on the 
site. The mutator mice lacking the 3 '-5' exonuclease accumulate 
mutations at a frequency that is in range of the mtRNA 
polymerase misincorporation. Gonsequently, we only obser\'e 
positive correlations between mtDNA mutagenesis and mtRNA 
integrity in these mutator mice. 

Our results are in line with deep-sequencing analysis of liver 
mtDNA [18] that questions the general effect of aging on mtDNA 
mutagenesis. As we have focused on seven 4-base sequences, we 
cannot rule out that mutations outside these sites contribute to the 
age-associated mitochondrial dysfunction. Our data nevertheless 
indicate that some sites are more prone to mutagenesis with age, 
and may correlate with those that undergo progressive clonal 
expansion to preclude mitochondrial dysfunction [2,4,19]. Our 
method can be used to identify the potential of a position to 
undergo clonal expansion. Garriers of the common deletion 
mutation may have up to 60% of mutated mtDNA molecules 
without manifesting into disease [12]. It would be interesting to 
address mtRNA integrity in carriers of the common deletion, to 
correlate this with the onset of disease, known as the Kearns-Sayre 
syndrome. 

In this study, we focused on brain mtDNA mutagenesis. Brain 
and heart were previously shown to accumulate mitochondrial 
mutations to a similar extent with age [9] . We previously detected 
equal mutation frequencies in the same sites in liver and brain, but 
higher in lung [15]. As a rough estimate, the average brain 
mtDNA mutation frequency in all seven sites is 0.7 x 10'* per nt, 
which is in coherence with results obtained by the deep- 
sequencing analysis of liver mtDNA [18]. Recent ultra-sensitive 
sequencing of human brain mitochondria has shown that age 
indeed increase mutagenesis. The apparent contradiction to our 
data could be that the relative age of the mice in our cohort was 
relatively much lower than the human cohort in the cited 
reference [20]. In contrast to mouse brain, the frequency of 
mtDNA mutations in liver did not change with age in any sites and 



suggests that age influences mtDNA mutagenesis in a tissue/cell 
type-specific manner. In particular, this is demonstrated by the 
reduction in mtDNA mutation level in colorectal cancer [21]. 

Interesting age-independent relationships between the mutation 
frequencies within individuals were discovered. We found a 
positive correlation between Ndl and Nd3, which might be 
expected as both sites showed age-associated mutagenesis. 
However, a stronger correlation was observed between mutation 
frequencies in 12S and Ndl, as well as in 12S and Nd6, which was 
independent of age. In comparison, none of the other sites (Nd5, 
Gytb, Coxl) seemed to be interconnected. Perhaps a subset of 
mtDNA genes/ sequences is particularly important in subregions of 
the brains, and that the quality of these regions is selected for. An 
analogy to this is observed during the purifying deselection of non- 
synonymous mutations in the protein-coding mtDNA genes 
during transmission [6], suggesting that mtDNA mutagenesis is 
under functional control. An analogue indication of a targeted 
mutagenesis was suggested in very young human brain samples as 
well [20]. 

The presence of nuclear mtDNA fragments (numts) [16] could 
theoretically introduce errors as we analyze total DNA and 
pseudogenes with altered nucleotide sequence potentially could 
generate positive signals in our restriction cleavage inhibition- 
based assay. However, since the primers used in the qPCR display 
similar efficiency and the GT values are comparable for the seven 
sites, it means that the mtDNA copy number is identical 
independent of the site to quantify. Again, this demonstrates that 
numts do not contribute significandy to the readout in any of the 
seven sites. Additionally, the melt-curve analyses do not indicate 
the presence of another amplicon to be amplied. Nevertheless, 
only the numts that deviate from mtDNA sequences would be 
detected, thus the potential outcome of numts on our conclusion 
would then be that the real mtDNA mutation rates are lower than 
we report. The lack of correlation between mutation frequency 
and mtRNA error frequency makes it very unlikely that numts are 
sufficiently frequent to interfere with our results. 

The expression of mitochondrial genes in the old cohort was not 
reduced, supporting our findings that mtRNA integrity was 
unaffected by aging as substitutions and deletions are likely to 
impair the stability of mtRNA. In view of our data we therefore 
believe that the correlation between mtDNA deletion and 
mitochondrial dysfunction with age is merely a secondary effect 
with minor functional impacts. In fact, the general expression of 
mitochondrial transcripts was slightly increased in old mice, 
presumably as an adaptive response to the observed mitochondrial 
function impairment. The innocuous impact of mtDNA mutations 
on age-related mitochondrial dysfunction is supported by the 
parallel reduction in complex II, which is entirely nuclear encoded 
(Figure SI). 

In summary, addressing mtRNA error frequency to evaluate 
mtDNA mutagenesis is a vital approach to determine the impact 
of mitochondrial mutations and can be used to test putative 
contribution of mtDNA mutations in various diseases. 

Supporting information 

File SI This Supporting Information file combines 
Supporting Methods and Table SI. Figure SI and S2 are 
provided as separate files. 
(DOGX) 

Figure SI Mitochondrial dysfunction with age. Enzymat- 
ic activity measurements of complex I, II and V in isolated brain 
mitochondria reveal an average 35% reduction in old (18 months. 



PLOS ONE I www.plosone.org 



6 



May 2014 I Volume 9 | Issue 5 | e96940 



mtRNA Errors Override mtDNA Mutations 



n = 11) compared to young mice (1 month, n = 4). Figure shows 
mean with SD, p**< 0.01. 

(PDF) 

Figure S2 Overestimate of 28S rRNA error frequency. 

Errors in nuclear 28S ribosomal RNA were estimated with a 
shghtly modified procedure. RNA was reverse transcribed, and the 
cDNA/RNA hybrid was treated with 100 units of TaqI for 1 h at 
65°C and followed by qPCR analysis. Figure shows mean with 
SD. 
(PDF) 



Acknowledgments 

We thank Alexander Rowe for critical reading of the manuscript. 

Author Contributions 

Conceived and designed the experiments: W^V KS YE JMS LE. 
Performed the experiments: WW KS YE JMS. Analyzed the data: WW 
KS YE JMS JBS MB LE. Contributed reagents/ materials/ analysis tools: 
JBS. Wrote the paper: WW KS JMS MB LE. 



References 

1. Linnane AW, Marzuki S, Ozawa T, Tanaka T (1989) Mitochondrial DNA 
mutations as an important contributor to ageing and degenerative diseases. 

Lancet: 642-645. 

2. Payne BA, Wilson IJ, Hatdcy CA, Horvath R, Santibanez-Korcf M, ct al (201 1) 
Mitochondrial aging is accelerated by anti-retroviral therapy through the clonal 
expansion of mtDNA mutations. Nat. Genet 43: 806-810. 

3. WanagatJ, Cao Z, Pathare P, Aiken JM (2001) Mitochondrial DNA deletion 
mutations eolocalize with segmental electron transport system abnormalities, 
muscle liber atrophy, fiber splitting, and oxidative damage in sarcopenia. 
FASEBJ 15: 322-332. 

4. CampbeU GR, Rraytsberg Y, Krishnan KJ, Ohno N, Ziabreva I, et al (2012) 
Glonally expanded mitochondrial DNA deletions within the choroid plexus in 
multiple sclerosis. Acta Ncuropathol 124: 209-220. 

,"). Wai T, Tcoh D, Shoubridge EA (2008) The mitochondrial DNA genetic 
bottleneck results fi-om replication of a subpopulation of genomes. Nat. Genet. 
40: 1484-1488. 

6. Stewart JB, Preyer C, Elson JL, Wredcnberg A, Cansu Z, et al (2008), Strong 
purifying selection in transmission of mammalian mitochondrial DNA. PLoS. 
Biol6; elO. 

7. Sproule DM, Kaufmann P (2008) Mitochondrial encephalopathy, lactic acidosis, 
and strokelike episodes: basic concepts, clinical phenotype, and therapeutic 
management of MELAS syndrome. Ann. N. Y. Acad. Sci. 1142: 133-158. 

8. Trifunovic A, Wredcnberg A, Falkenberg M, Spelbrink JN, Rovio AT, et al 
(2004) Premature ageing in mice expressing defective mitochondrial DNA 
polymera.sc. Namre 429: 417 -423. 

9. Vermulsl M. BirlasJH. Kujolh CJC, Lariigrs WG. Raljinoxilch PS. c\ al (2007). 
Mitochondrial point mutations do not limit the natural lilespan of mice. Nat. 
Genet. 39: 540-,S43. 

10. Vermulst M, WanagatJ, Kujoth GG, Bielas JH, Rabinovitch PS, et al (2008) 
DNA deletions and clonal mutations drive premature aging in mitochondrial 
mutator mice. Nat. Genet. 40: 392-394. 



11. McDonald DG, McMcnamin JB, Farrell .MA, Droogan O, Green AJ (2002) 
Familial childhood onset neuropathy and cirrhosis with the 4977bp mitochon- 
drial DNA deletion. .Am. J. Med. Genet. Ill: 191-194. 

12. Nakada K, Inotie K, Ono T, Isobe K, Ogura .\, et al (2001) Inter-mitochondrial 
complementation: Mitochondria-specilic system preventing mice from expres- 
sion of disease phenotypes by mutant mtDNA. Nat. Med. 7: 934—940. 

13. Edgar D, Shabalina I, Camara Y, Wredcnberg A, Calvaraso MA, et al (2009) 
Random point mutations with major effects on protein-coding genes are the 
driving force behind premature aging in mtDNA mutator mice. Cell Metab 10: 
131-138. 

14. Kraytshcrg Y, Simon DK, Turnbull D.\I, Khrapko (2009) Do mtDNA deletions 
drive premature aging in mtDN.A mutator mice? Aging Gell 8: ,502-506. 

15. Halsne R, Esbensen Y, Wang W, Schdiler K, Suganthan R, et al (2012) Lack of 
the DNA glycosylases MYH and OGGl in the cancer prone double mutant 
mouse does not increase mitochondrial DNA mutagenesis. DNA Repair (Amst) 
11: 278-285. 

16. Ricchetti M, Tekaia F, Dujon B (2004) Gontinued colonization of the human 
genome by mitochondrial DNA. PLoS. Biol. 2: E273. 

17. Kuhn CD, Geigcr SR, Baumli S, Gartmann M, Gerber J, et al (2007) Functional 
architecture of RNA polymerase 1. Gell 131: 1260-1272. 

18. Ameur A, Stewart JB, Freyer C, Hagstrom E, Ingman M, et al (201 1) Ultra-deep 
sequencing of mouse mitochondrial DNA: mutational patterns and their origins. 
PLoS. Genet. 7: el002028. 

19. Ishikawa K, Takenaga K, Akimoto M, Koshikawa N, Yamaguchi A, et al (2008) 
ROS-generating mitochondrial DNA mutations can regulate tumor cell 
metastasis. Science 320: 661—664. 

20. Kennedy SR, Salk JJ, Schmitt MW, Loeb LA (20 1 3) Ultra-sensitive sequencing 
reveals an age-related increase in somatic mitochondrial mutations that are 
inconsistent with oxidative damage. PLoS. (ienet. 9: el003794. 

21. Ericson NG, Kulawiec M, Vermulst M, Sheahan K, O'SuUivan J, et al (2012) 
Decreased mitochondrial DNA mutagenesis in human colorectal cancer. PLoS. 
Genet. 8:cl002689. 



PLOS ONE I www.plosone.org 



7 



May 2014 I Volume 9 | Issue 5 | e96940 



